Quasi-energies and Floquet states of two weakly coupled Bose-Einstein condensates 

under periodic driving 
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We investigate the quasi-energies and Floquet states of two weakly coupled Bose-Einstein conden- 
sates driven by a periodic force. The quasi-energies and Floquet states of this system are computed 
within two different theoretical frameworks: the mean-field model and the second-quantized model. 
The mean-field approach reveals a triangular structure in the quasi-energy band. Our analysis of 
the corresponding Floquet states shows that this triangle signals the onset of a localization phe- 
nomenon, which can be regarded as a generalization of the well-known phenomenon called coherent 
destruction of tunneling. With the second quantized model, we find also a triangular structure in 
the quantum quasi-energy band, which is enveloped by the mean-field triangle. The close relation 
between these two sets of quasi-energies is further explored by a semi-classical method. With a 
Sommerfeld rule generalized to time-dependent systems, the quantum quasi-energies are computed 
by quantizing semiclassically the mean-field model and they are found to agree very well with the 
results obtained directly with the second-quantized model. 

PACS numbers: 03.75.Kk, 03.75.Lm, 42.50.Vk, 05.30.Jp 



I. INTRODUCTION 



Due to its simplicity, a single particle in a double- 
well potential has been a paradigm to demonstrate many 
fundamental quantum phenomena, in particular, quan- 
tum tunneling and its control[l|. Immediately after 
the experimental creation of Bose-Einstein condensates 
(BECs) with dilute alkali atomic gases people re- 

alize the new possibility of putting a BEC in a double- 
well potential and using it to mimic this paradigm sys- 
tem to demonstrate experimentally quantum tunneling 
and other fundamental quantum phenomena. The sub- 
sequent studies show that a BEC in a double- well poten- 
tial has richer physics due to interaction. For example, it 
was found that the tunneling of BEC between the wells 
can be suppressed and therefore self-trapped in one of the 
wells[3|-[5|. This self-trapping phenomenon has now been 
observed experimentally with a BEC[1,Q- More interest- 
ingly, the nonlinear two-mode model derived to describe 
a BEC in a double-well potential was found to be able to 
describe the tunneling between Bloch bands for a BEC 
in an optical lattice 8]. Due to interaction, a new quan- 
tum phenomenon called nonlinear Landau-Zener tunnel- 
ing was predicted and later observed in experiment [^Q. 

It is known that, for a single particle in a double- well 
potential, one can use an external periodically driving 
field to control quantum tunneling, either enhancing [lo|- 
or suppressing it[l^-[2^ One then wonders whether 
this kind of control can be also achieved for a BEC in 
a double-well potential. There have been several stud- 
ies in this regardjlJ-flOl- These studies indeed find that 
the periodically driving force can strongly affect the tun- 
neling between two weakly coupled BECs and therefore 
be used to control the tunneling. Recently, we found 
that such a control of quantum tunneling can also be 
achieved in an optical waveguide system^ and be used 
to improve the performance of an all-optical switch [33|. 



In this paper we investigate the quasi-energies and Flo- 
quet states of two weakly coupled BECs under periodic 
driving, which can be realized experimentally with either 
a double- well potential or an optical lattice [3^. Quasi- 
energies and Floquet states are two basic concepts and 
tools in describing and understanding periodically driv- 
ing systems. One can use either a mean-field nonlinear 
two-mode model or a second quantized model to describe 
such a system. In this paper we use both models to 
compute the quasi-energies and Floquet states. In the 
mean-field two-mode model, we discover that there can 
be more than two Floquet states and quasi-energies in a 
certain range of parameters that characterize the ampli- 
tude and frequency of the modulating force. With these 
additional Floquet states, there appears a triangle in the 
quasi-energy levels. This triangular structure in quasi- 
energies turns out to be crucial to understanding the lo- 
calization phenomenon that has been found and studied 
previously [24, 25, 27]. Our analysis shows that the local- 
ization phenomenon can be regarded as a generalization 
of a well-known phenomenon called coherent destruction 
of tunneling(CDT). Therefore, we call it nonlinear coher- 
ent destruction of tunneling(NCDT)[3l|. 

In the second quantized model, our computation also 
reveals a triangular structure in the quasi-energy levels. 
Interestingly, the quantum triangle is enveloped perfectly 
by the mean-field triangle, indicating a close connection 
between these two different approaches. By analyzing 
the corresponding Floquet states, we find that this quan- 
tum triangle of quasi-energies is also connected to the lo- 
calization phenomenon called NCDT. The close relation 
between quantum quasi-energies and mean-field quasi- 
energies is further explored by a semi-classical method. 
By using a Sommerfeld quantization rule adapted for 
a time-dependent system, we re-calculate the quantum 
quasi-energies by quantizing semiclassically the mean- 
field model. The results match very well with the quan- 
tum quasi-energies obtained by directly using the second 
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quantized model. 

Due to the complication brought by the chaos in the 
region of moderate frequencies, the focus of our paper is 
on cases of high frequency modulation. 



II. QUASI-ENERGIES AND FLOQUET STATES 

We consider a system of N identical bosons, which can 
occupy only two quantum states. If there is interaction 
between bosons, the system Hamiltonian reads [1] 



2N 



{a^'a^aa + b^b^bb) , 



(1) 



where 7 is the energy difference between the two quan- 
tum states denoted by a\a and w,b and v is the cou- 
pling constant between the two modes. The interaction 
strength is given by 



(2) 



where we have used a reasonable assumption that the 
wave functions of the two quantum states are the same 
except a possible trivial shift of the center and the wave 
function is normalized J \ il;Q{r)\'^d'r = 1. 

When the temperature is very low so that we can ig- 
nore any thermal effect and at the same time the number 
of bosons N is very large, it is appropriate to make the 
following coherent substitutes 



(3) 



a={a)/VN, b={b)/VN. 



This leads to a mean-field Hamiltonian 



N 



(4) 



The system described above has now been realized with 
a double- well potential. For the experiment in Ref.[l], 
there are about 1150 atoms and a simple estimate gives 
V « 65.3s~^ and c/v k, 15. This system can also be 
realized experimentally with an optical lattice [sl. |33|. 

In our study, we have 7 = A cos(a;t), that is, the energy 
difference between the two quantum states is changed 
periodically. With the double-well potential, this can be 
achieved by shifting periodically the power of lasers that 
generate the double wells. For an optical lattice, this can 
be accomplished by shaking along the lattice direction. 
We focus our study on the quasi-energies and Floquet 
states of this system as these are two basic concepts and 
tools in understanding a periodically driving system. 



A. Mean-field model 

We first consider the mean-field model. From the 
mean- field Hamiltonian ([?]), we can obtain a two- mode 
Gross-Pitaevskii equation 



dt 
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(5) 



where we have used the natural unit h — 1. Although the 
parameters c, w, and oj are of unit of energy, we shall 
treat them as dimensionless parameters in the following 
discussion because what is essential is the ratios between 
these parameters not their absolute values. 

Like its linear counterpart, a nonlinear periodic time 
dependent equation admits solutions in the form of Flo- 
quet states. For Eq.®, its Floquet state has the follow- 
ing form 



-3 — Z£t 



(6) 



where both (t) and (j)2 (t) are periodic functions of pe- 
riod of T = 2tt /lo and the constant e is the corresponding 
quasi-energy. After one period, this solution returns to 
its original state by picking up an extra phase of eT. To 
calculate numerically Floquet states and quasi-energies, 
we follow the strategy that was used to compute non- 
linear Bloch states and the eigen-energies (sif . In this 
strategy, we expand the Floquet states in Fourier series 



n— — L 



71 — — L 



(7) 



where L is the cut-off and equal to 10 in our computa- 
tion. With the substitution of the above Fourier series 
into Eq. (O, one can obtain AL + 2 equalities for the co- 
efficients of each Fourier term e*"'^*. The Floquet state 
and the quasi-energy are found by finding the roots of this 
set of 4L + 2 nonlinear equations. Our method is differ- 
ent from the previous methods used to compute Floquet 
states and quasi-energies [24j . We believe that it is more 
powerful. For example, it can find the Floquet states that 
correspond to hyperbolic fixed points in Poincare section, 
which can not be found with the previous method [2^. 

Our numerical results of quasi-energies are plotted in 
FigU It is clear from FiglT] that, for the linear case, 
there are two quasi-energies at a given value oi A/ui with 
one isolated degeneracy point. For the nonlinear case, we 
notice that there are three quasi-energies within a certain 
range of A/uj with two of them degenerate. The three 
quasi-energies form a triangle in the quasi-energy levels 
as seen in Figmb). Among the three quasi-energies, two 
quasi-energy levels are similar to their linear counterparts 
with one isolated degenerate point while the third quasi- 
energy level has no linear counterpart. Moreover, the 
third quasi-energy is degenerate and corresponds to two 
different Floquet states; this is indicated by marking the 
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same point in Figlljb) with two symbols Pi and P2 ■ Note 
two things: (1) there is no threshold value of c for the 
triangle to appear; (2) the right corner of triangle is open 
for relatively larger nonlinear parameter c. 
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FIG. 1: (Color online) Quasi-energies as a function of Ajio at 
(a)c = 0; (b)c = 0.4; (c)c = 0.8. Solid lines are numerical 
results and circles are the approximate analytical results for 
high frequencies with Eqs. (|9ll0|) . v = = 

Despite the obvious similarity between the nonlinear 
Floquet states and the linear ones, there are a couple 
of conceptual differences. (1) A periodically driven n- 
level linear system possesses precisely n Floquet states 
whereas the number of nonlinear Floquet states of n- 
mode system can be bigger than n as we have witnessed 
above. (2) In linear case, all wave functions can be de- 
composed into a superposition of Floquet states and. 



therefore, the dynamics of the system is dictated by Flo- 
quet states. In the nonlinear case, the superposition prin- 
ciple breaks down, the dynamics of the system can no 
longer be completely determined by Floquet states. 

The triangular structure of the quasi-energy is very 
similar to the energy loop discovered within the context 
of nonlinear Landau-Zener tunneling [§!]. In fact, they 
are mathematically related. For high frequencies, w ^ 
max{w, c}, we take advantage of the transformation 



a exp[ 



-I ], 6 = 5exp[i ^ ]. (8) 



2tj 



After averaging out the high frequency terms l2£ 
obtain a non-driving nonlinear model. 



^ijo(AMfe' + c|a'|V, 



ih' 



-Jo(A/tj)a' + c|6'|26' 



IJ, we 

(9) 
(10) 



where Jq is the zeroth-order Bessel function. It is clear 
from the transformation in Eq.([5]) that the eigenstates 
of the above non-driving nonlinear equations correspond 
to the Floquet states of Eq.([5]). We have computed the 
eigenstates of Eas. (|9ll0p and the corresponding eigenen- 
ergies, which are plotted as circles in Fig[T] The consis- 
tency with our previous numerical results is obvious. As 
is known in Ref . [8] , the above nonlinear model admits ad- 
ditional eigenstates when c > Jq{AIijj)v. Therefore, this 
can be regarded as the condition for the extra Floquet 
states to appear for the driving nonlinear model Eql5] at 
high frequencies. Since the Bessel function Jo(A/w) can 
be zero, there is no threshold value of c for the triangle 
to appear in the quasi-energy band. 

The nonlinear Floquet states are also examined thor- 
oughly. We find that some of them are localized, which is 
very different from the linear Floquet states that are al- 
ways unlocalized. To describe localization, we introduce 
a new variable, p = (|ap — |5p)/2, which measures the 
population difference between the two modes. One Flo- 
quet state is localized if the average of p over one period. 



1 

T 



(11) 



is nonzero; it is unlocalized if (p)t = 0. In Figl2] , the 
population difference p is plotted as a function of time for 
three stable nonlinear Floquet states marked as Pi , P2 , P3 
in Figlljb) . Evidently, one of these states is unlocalized 
since p oscillates around zero. However, two other states 
are localized with p oscillating around a non-zero value. 
The localization means that the BEG described by such 
Floquet states tends to stay in one mode and reluctant 
to tunnel to the other mode. Therefore, localization can 
be understood as a suppression of tunneling. Our study 
shows that on one hand, all the localized Floquet states 
correspond to the highest quasi-energies on the triangle 
and on the other hand, all Floquet states in the linear 
case and all the Floquet states not related to the quasi- 
energy triangle are not localized. This implies that the 
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triangle in Figlljis related to localization or suppression structure in quasi-energies are related to localization or 

of tunneling. This is indeed the case as we have shown suppression of tunneling, 
in Ref-fsij. We shall not repeat what we have done in 
Ref. [3l| ; we shall look into this connection from a differ- 




FIG. 2: Population imbalance p for three stable nonlinear 
Floquet states at c = 0.4,?; = = IQ, A/lj = 2.0, for an 
interval of four periods of the driving force(solid lines). The 
Floquet states correspond to the quasi-energies in Fig[llb) 
marked as Pi,P2,Pz with triangles. The squares are for Flo- 
quet states in the highest two quantum quasi-energy levels 
with iV = 500. 

In the mean-field model ([51), the norm |ap -|- |6p = 
1 is conserved and the overall phase is not essential to 
the dynamics. Therefore, we can reduce the complex 
dynamical variables a = loje*^", b = \b\e^^'' to a pair of 



real variable, p — 



-|6p)/2 and the relative phase q — 



9b — Oa- In terms of p and q, the mean-field Hamiltonian 
dH) becomes 



Hci = Apcos{Ljt) - -Vl~V cosq+ -(4p2 + 1). (12) 

As p and q are canonically conjugate variables of the 
above classical Hamiltonian system, one can derive a set 
of equations of motion. From the equations of motion, 
one can plot the Poincare section of this system. Two 
Poincare sections are illustrated in Fig[3] for two sets of 
parameters. As the overall phase is removed, the Floquet 
states correspond to the fixed points in Poincare section. 

The parameters for FiglSl^a) are outside the triangle 
range. In this figure, there are only two fixed points 
located at p = and all the motions around the fixed 
points are oscillating around p — 0, indicating no lo- 
calization or suppression of tunneling. The situation is 
different in Fig[3I^b), whose parameters lie in the triangle 
range. In Figl3{b), there are four fixed points: one at 
q = 0(or 27t); three at g = tt. Among the three at g = tt, 
one is hyperbolic and unstable whereas the other two are 
not only stable but localized. Moreover, all the orbits 
surrounding these two stable fixed points at q = n are 
localized solutions. These again show that the triangle 



In Ref.[3l|, the localization phenomenon discussed 
above is called nonlinear coherent destruction of tun- 
neling (NCDT). There are two reasons for this. First, 
the degeneracy point in FiglIJa) is related to a local- 
ization phenomenon called coherent destruction of tun- 
neling (CDT) and the triangle can be seen as the re- 
sult of enlargement of the degenerate point by nonlin- 
earity. Second, as we have seen in FiglJl the localiza- 
tion phenomenon is intimately related to the nonlin- 
ear Floquet states and we know that CDT is related 
to linear Floquet states. The localization phenomenon 
which we call NCDT has been called in literature self- 
trapping or, more precisely, periodically modulated self- 
trapping [11 [H, [27 1 . 



B. Second quantized model 

We now turn to the second quantized model H]) and 
compute its Floquet states and quasi-energies. For a non- 
driving system, it is well known that the eigen-energies 
and eigenstates of the second quantized model are closely 
connected to its mean- field counterparts [sgI Ist} . For this 
periodically driving system, we want to explore how its 
quantum Floquet states and quasi-energies are related 
to its mean-field counterparts and the localization phe- 
nomenon called NCDT. 

We follow the well established Floquet theory for a 
quantum system [ssl [39j to compute numerically quan- 
tum Floquet states and quasi-energies. In the pro- 
cess, we have converted the second quantized Hamilto- 
nian (HI) into a pseudo-spin Hamiltonian by introducing 
three angular momentum operators Jx = (a^6 + b^a)/2, 
Jy = i{b^a — ba^)/2, and Jz = {a)a — ¥b)/2, for which the 
Casimir invariant is = {N/2){N/2 + 1). The second 
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FIG. 4; Quantum quasi-energies (A'^ — 40) as a function of 
A/lo at V = 1,uj ^ W for (a) c = 0.0; (b)c = 0.4; (c)c = 
0.8. The open circles are mean-field quasi-energies. Note that 
for comparison with mean-field theory, the quantum quasi- 
energies have been divided by A''. 



spanned by TV -f 1 spin states | J = N/2, = M) with 
M = -N/2, -N/2 N/2. 

Our numerical results for quantum quasi-energies for 
= 40 are shown in FigUl We immediately notice 
that these quantum quasi-energy levels have very sim- 
ilar structures to their mean-field counterparts. For the 
non-interacting case in FigHJa), there is a single degen- 
eracy point. For interacting cases in Fig|31[b)&(c), there 
are triangular structures just as in the mean-field model. 
For comparison, the mean-field quasi-energies are plotted 
as open circles in FigU) To one's amazement or expec- 
tation, the quantum quasi-energies are bounded by the 
mean-field results perfectly. Another interesting feature 
in Figgis that all the quasi-energies in the triangle area 
is doubly degenerate and this degeneracy immediately 
breaks up outside the triangle. The feature is related the 
localization phenomenon NCDT as we shall discuss next. 

There is also a close relation between quantum Flo- 
quet states and mean-field Floquet states. We examine 
this relation in terms of localization. To measure how a 
quantum Floquet state is localized, we define 



{P)t 



1 

T 



dt{Ur^{t)\J,\Ur^{t)) 



(14) 



for a given Floquet state |u„(t)). This variable {P)t 
quantifies the population difference between the two 
modes. We have plotted this variable for certain quan- 
tum Floquet states in Fig. O It is apparent from this 
figure that only the Floquet states for the quasi-energies 
inside the triangle are localized. This again establishes 
the connection of the triangle (quantum or mean-field) 
to the localization phenomenon NCDT. This localization 
also explains why the Floquet states inside the triangle 
are doubly degenerate. When localization occurs, there 
are two equal possibilities. It can localize either in mode 
a or in mode b; this leads to degeneracy. The mean-field 
results are also plotted in Fig. O They match very well 
with the results for the two highest quantum Floquet 
states. This good correspondence can be more clearly 
seen in Fig. [21 the temporal evolution of two highest 
quantum Floquet sates agrees very well with the mean- 
field results for an interval of four periods of the driving. 

The quantum quasi-energies and Floquet state were 
studied in Ref. [l^l ■ Their relation to the localization was 
also examined there. Our primary purpose here is to 
compare them to the mean-field results and explore their 
relations, which has not been studied so far. 



quantized Hamilton of the system then becomes 

Hq = -vJ^ + if + Acos{ujt)J,_ + ^{N - 2). (13) 

With this transformation, our system of N identical 
bosons becomes a spin system, whose Hilbert space is 



III. SEMICLASSICAL QUANTIZATION 

In the previous section, we have demonstrated by di- 
rect numerical computation how the quantum Floquet 
states and quasi-energies are connected to their mean- 
field counterparts. This relation can be further explored 
with a semiclassical method as the mean-field model ^ 
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FIG. 5: Population difference {P)t for every Floquet state 
in tlie fiigliest two quantum quasi-energy levels(solid line), 
the 149th and 150th quantum quasi-energy levels(dot-dashed 
line), and the 249th and 250th quantum quasi-energy lev- 
els(dotted line) at c/v = 0.4, cu/v = W,N = 500. The open 
circles are for the population difference {p)t for the highest 
mean-held quasi-energy level in Figlljb). 

can be regarded as the classical limit of the second quan- 
tized model ([IJ in the limit of TV ^ (X) (4 C | . We shall 
follow the procedure in Ref. [13, Sll, l42l |43|| and try 
to quantize the classical Hamiltonian in Eg. p^ . which 
is equivalent to Hamiltonian ([4]), with the Sommerfeld 
rule. However, as our system is time dependent, the usual 
Sommerfeld quantization rule has to be generalized. 




FIG. 6: Periodic vortex tube. Two paths are shown. The 
path 7i lies in a plane of f=const and 72 is a path connecting 
a point (p, q) at time t with the same point at time t + T. 

The generalization of the Sommerfeld rule has been 
done for any time-dependent systemjH, |4^. The basic 
idea is to regard time as a dynamic variable and intro- 
duce a new canonical momentum which conjugates time. 
We shall not go into the details of this theory and shall 
only describe how this generalization works for the case of 
our interest, a periodic time dependent system. As seen 
in the Poincare section of FigO there are closed orbits 
around fixed points. These closed orbits will change their 
positions and shapes in the phase space with time and re- 
turn to their original points and shapes after one period. 
This kind of evolution forms a tube in the space spanned 



by p, q, t as depicted in FiglB] This tube is called vor- 
tex tube. As the system is periodic in time, the tube in 
FiglB]is essentially a torus. The quantization can be done 
by choosing two independent closed paths on the vortex 
tube which cannot be homotopically deformed onto each 
other and requiring 

/i = pdq^nih/N, (15) 

If T 

h = TT f (pd'Q - Hcidt) + —e 

= n2h, (16) 

where ni and 712 non-negative integers. The quantization 
is done in two steps: (1) we first find a path 71 that fulfills 
the quantization condition for /i; (2) the quantization 
condition for I2 is then used to compute the quasi-energy 
e as 

eni.na = - f {pdq - Hddt) + n2U} . (17) 

In the above, n2Lo means that quasi-energy e is only de- 
fined modulo w, refiecting the unique nature of quasi- 
energy. One can view h/N in Eq. fTS)) as the effective 
Plank constant [U nil, [4^, which goes to zero at the limit 
of iV ^ 00. 

Our semiclassical results of quasi-energies are plotted 
in FigOto compare with the quantum quasi-energies ob- 
tained directly from the second quantized model. They 
match perfectly, indicating the success of the generalized 
Sommerfeld quantization rule. In our calculation, the 
path 71 is chosen as the closed orbit in the Poincare sec- 
tion and 72 is the path along the maximal points of p on 
the tube as illustrated in FiglHl Note that the natural 
unit fi = 1 is used in our calculation. 

These semiclassical results are very helpful in under- 
standing why the quantum quasi-energies are enveloped 
by the mean-field quasi-energies as seen in FigHl We first 
look at the simple case when there are only two fixed 
points in the Poincare section, as in Figl3fa). The fixed 
point aX q — Q corresponds to the nonlinear Floquet state 
with lower quasi-energy and the other fixed point cor- 
responds to the Floquet state with higher quasi-energy. 
This implies that the quantization for orbits around the 
fixed point at g = produces quasi-energies that are 
higher than the corresponding mean-field quasi-energy 
and the quantization for orbits around the fixed point 
at g = TT yields quasi-energies that are lower than the 
corresponding mean-field quasi-energy. As a result, the 
quantum quasi-energies are bounded by the mean-field 
quasi-energies. The double degeneracy of the quantum 
quasi-energies within the triangle can also be explained 
with this semiclassical approach. As shown in Fig[3Kb), 
there are two stable fixed points aX q — it. These two 
fixed points correspond to two Floquet states with the 
same quasi-energy. This indicates that if one quantizes 
semiclassically the orbits around these two fixed points, 
one would get two identical sets of quasi-energies. This 
explains the double degeneracy. 
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FIG. 7: Comparison between the quantum quasi-energy levels 
(solid lines) with iV = 40 and semiclassical quasi-energy levels 
(open circles) at c = 0.4, w = = 10. (a) nondcgcncratc 
quasi-energy levels; (b) degenerate quasi-energy levels. For 
clarity, we have only plotted a portion of the quasi-energy 
levels. 



IV. CONCLUSIONS 



To summarize, wc have stiidied the quasi-energies and 
Floquet states of two weakly coupled Bose-Einstein con- 
densates subject to a periodic driving. Both the mean- 
field model and the second quantized model are used. A 
triangular structure was found in both mean-field quasi- 
energy levels and quantum quasi-energy levels. More- 
over, we have revealed that the quantum quasi-energy 
levels are bound by their mean-field counterparts and we 
have explained it with semiclassical quantization. In ad- 
dition, by looking into the Floquet states, we have found 
that the triangle in the quasi-energies is related a lo- 
calization phenomenon which we call nonlinear coherent 
destruction of tunneling (NCDT). 
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